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The scaling of the viscosity of polymer melts is investigated with regard to the molecular weight. 
We present a generalization of the Rubinstein-Duke model, which takes constraint releases into 
account and calculate the effects on the viscosity by the use of the Density Matrix Renormalization 
Group (DMRG) algorithm. Using input from Rouse theory the rates for the constraint releases are 
determined in a self consistent way. We conclude that shape fluctuations of the tube caused by 
constraint release are not a likely candidate for improving Doi's crossover theory for the scaling of 
the polymer viscosity. 



I. INTRODUCTION 

The visccjsity of polymer melts has been investigated 
intensively.!!! The models for describing this behavior, 
however, are not yet satisfying. Experiments show that 
the viscosity rj scales like where M is the 

molecular weight. This behavior is valid for several 
decades of the molecular weight. An early approach 
was the reptation model by dc GenncM which yields 
r] oc in the lirait of infinitely long polymers. For 
short polymers Doic! calculated the viscosity by taking 
tube length fluctuations into account and found a region 
of M where the scaling is of the correct size - but this 
region appears to be much too small. A discrete model 
for reptation which includes tube length fluctuations is 
the Rubinstein-Duke (RD) model.!3'E! Recently it was 
shownQ that this model does not only provide a region 
where a scaling of iVI? '^*°-^ can be found, as already 
shown by RubinsteinJ^ but also that this model shows 
a crossover to the reptation exponent 3. But again the 
mass region where the correct exponent is valid does not 
exceed one order of magnitude. In this paper we extend 
the RD model in order to investigate whether constraint 
release (CR) broadens that region. While good models 
for the non-linear regime of viscosity can be constructed 
by using CR,0'E! there was only little success in the linear 
regimes 



II. THEORY 

A. RD model 

In order to investigate the role of tube length fluc- 
tuations in reptation theory Rubinstein introduced a 
discrete "repton" model which allows the description 
of three-dimensional reptation dynamics by a one- 
dimensional lattice gas.Q Duke generalized the model 
to the case thai, an external electric field acts on a 
charged polymem which allows a good description of thft 
diffusion constant in gel electrophorese experimcnts.t2l 
However, also in the absence of a field the generalization 
by Duke is useful in so far as it provides a reference axis 



along which the displacement of the polymer chain as 
a whole can be monitored. This makes it possible to 
calculate the diffusion coefficient and the viscosity from 
the model, without resorting to independent hypotheses. 
Rubinstein assumes that the constraints of the other 
polymers divide space into cells which form a d- 
dimensional regular cubic lattice. The polymer occupies 
a series of adjacent cells, the "primitive path". It is not 
possible for the polymer to traverse the edges of the cells 
(in two dimensions: the lattice points) so that only the 
ends of the polymer can enter new cells. The polymer is 
divided into segments whose length is of the order of the 
lattice constant, the number of segments is proportional 
to the length of the polymer or the molecular weight. 
The orientation of the lattice is introduced by Duke in a 
way that the electric field is diagonal to the lattice, i.e. 
in three dimensions in the (lll)-direction. 
A segment - called a "repton" - is allowed to jump into 
an adjacent cell according to the following rules: 

1. The reptons in the bulk are only allowed to jump 
along the primitive path. 

2. No cell in the interior of the primitive path may be 
left empty. 

3. The ends move freely insofar rule 2 is respected. If 
an end repton occupies the cell alone, it can only 
retract in the cell of the adjacent repton. If the 
adjacent repton is in the same cell, the end repton 
may enter any of the 2d surrounding cells. Rep- 
tons in the bulk jump with the same probability as 
reptons at the end into occupied cells. 

Rule 1 ensures that the polymer does not traverse the 
edges of the cells. Rule 2 is motivated by the fact that 
the segments are of the size of the lattice constant. Fi- 
nally rule 3 reflects the fact that there are more free ad- 
jacent cells for an end repton than occupied ones. If one 
considers a field F ^ the ratio of the probabilities to 
jump in respectively against the direction of the field is 
proportional to the Boltzmann weight. In the following 
we will only focus on the case without field, we only need 
the reference axis. 

As the shape of the primitive path is not affected by 
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FIG. 1: a) The repton model in two dimensions: the circles 
represent the reptons; the primitive path is marked by the 
bold lines, b) Projection onto one dimension. 



the movement of the polymer and only the ends are cre- 
ated or annihilated this model can be mapped onto one 
dimension. To this end a particle is assigned to each 
bond between reptons; starting on one end, a bond into 
the direction of the field is identified with an A-particle, 
against the field with a B-particle and finally a bond 
without change of potential, i.e. a bond between two 
reptons in one cell, by a vacancy (Fig. 
These particles are residing on a chain whose number of 
sites N is the number of bonds, thus the number of rep- 
tons minus one. As for each site three states are possible 
this rnodel can be identified with a "quantum" spin-one 
chain.BE3 We define 
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We label this set of states with X — {A, 0, B}. A state 
of a chain of length N can be considered as an element 
\ri) of the tensor base X — X®^; \r]) is constructed by 
the tensor product of the three component vectors for the 
individual sites. 

The dynamics of this one dimensional model is: 

1. A-particles may exchange with vacancies 0: A0 ^ 
0A. 

2. _B -particles may exchange with vacancies 0: 
B0 ^ 0B. 

3. A-particles may not exchange with i?-particles. 

By these rules it is guaranteed that the shape of the 
primitive path is conserved. The next rules define the 
boundary dynamics: 

4. At the ends of the chain A-particles may be anni- 
hilated: 

A-^ 0. 

5. At the ends of the chain B-particles may be anni- 
hilated: 

B^ 0. 



6. At the ends of the chain A-particles may be cre- 
ated: ^ A. 

7. At the ends of the chain i?-particles may be cre- 
ated: ^ B. 

The processes 1, 2, 4 and 5 take place with the same 
probability; the processes 6 and 7 are d-times more 
probable. 

In this way a stochastic interacting particle system 
on a one dimensional chain has been defined. As the 
transitions are independent of the previous history this 
process is Markovian. 



B. Quantum Hamiltonian 

A convenient way to describe the process mathemat- 
ically is the quantum Hamiltonian formalism which we 
will present here shortly, for details see Ref. HI 
The probability to be in state \ri) at time t is labeled by 
Pq{t). These probabilities of the individual states can 
be combined to a vector: \P{t)) = Due to 

the conservation of probability the entries of the vector 
\P{t)) sum up to 1 at any time t. With this definition 
the master equation can be written as: 

±\P[t))=-H\P{t)), (1) 

with the stochastic generator 

i? = - ^ ^ w,Y^^ \rj) (jy'l + XI XI 1*^^ ^'^l • 



(2) 

Here Wrj'^n is the transition probability from state \r]') 
to |ry). In other words the off diagonal elements of the 
matrix H are the negative transition rates between the 
respective states and the diagonal elements are the sum 
of the rates leading away from the respective state. 
The creation operators and are defined by a^0 = A 
and &^0 — B, acting on any other state yields zero. The 
annihilation operators a and b are defined by aA = 
and bB = 0, again acting on any other state yields again 
zero. Finally we define the number operators n"* = a^a, 
^ b^b and n'^ ^ 1 ~ — . By these definitions 
the stochastic generator H of the RD model reads:ll3 



N-l 

H = bi{d) + bN{d) + X 



(3) 



with 



6„(d) = d[n^ ~~al + n^ - b^] + ~ a„ + - 

Again, d labels the lattice dimension. We have chosen 
the time scale such that the hopping rate in the bulk 
equals unity. 
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C. Calculation of the viscosity 

The viscosity is proportional to the longest relaxation 
time of thep8tochastic generator H which is the inverse 
energy gap.tl Due to the conservation of probability the 
ground state of a stochastic generator has the eigenvalue 
zero, so that we only need to calculate the first eigen- 
value. 

A very efficient algorithm to calculate the lowest excita- 
tions of quantum spin chains is the densibt matrix renor- 
malization group (DMRG) algorithm.BEll For the stan- 
dard RD modd the usefulness of this algorithm has been 
dcmonstratedjj so that we employ this algorithm as well 
for the modifications of the RD model, to be introduced 
in what follows. 

The Hilbert space of the spin chain grows exponentially 
with the number sites TV; the number of states is 3^, 
so that for TV = 50 — the maximum number of sites 
considered in our calculations — the number of states 
is of the order 10^'^. By this it is reasonable to project 
all quantities onto a subspace consisting of the most im- 
portant states. The difficulty is to find out which states 
are the "most important". The DMRG algorithm is a 
method which provides a choice of states which is opti- 
mal in terms of a maximum of the probability with which 
the states contribute to the target state.EJ This proba- 
bility is gained by diagonalizing the density matrix. As 
the diagonalization of the density matrix for the whole 
system would be as laborious as the direct diagonaliza- 
tion of the hamiltonian the system is build up stepwise. 
Starting with a small system (e.g. 2 states) one adds 
iteratively two states until the system has reached the 
searched size. In each step the Hilbert space is reduced 
to the subspace of the most probable states, so that the 
system size increases while the dimension of the matrices 
remain constant. 



D. Constraint release 

The RD model as well as standard reptation theory is 
based on the assumption that the polymer moves in a 
fixed network formed by the surrounding polymers. This 
assumption is at best justifiable for a single polymer im- 
mersed in a gel. But for polymer melts, for which the 
viscosity is measured, it should be taken into account 
that the surrounding polymers reptate themselves. In 
terms of the Rubinstein model this means that the lat- 
tice itself is subjected to fluctuations. 
We consider the following model of lattice fluctuations: 
Imagine that a constraining polymer moves so far that 
the constraint for the investigated polymer is released so 
that it can move freely in this region. After a short time 
the constraining polymer returns or an other polymer has 
taken its place so that the free movement in this region 
is again prevented. The lattice has regained its originally 
structure but the primitive path may have changed in the 
bulk (Fig. 1^) . This is what we call a "constraint release" 





FIG. 2: A constraining polymer vacates for a short time its 
site and thus makes a constraint release event possible. 
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FIG. 3; Constraint release by ^_B-permutation in the RD 
model. 



(CR) event. 

A straightforward implementation of this mechanism is 
shown in Fig. ^. The turning of bend of the primitive 
path equals to the permutation of an AB-pair. At the 
ends of the path A-particles can be transformed into B- 
particles and vice versa which corresponds to a CR event 
at the ends. 

To construct the extended Hamiltonian which takes these 
mechanisms into account we introduce an operator which 
transforms B-particles into A-particles: 



— a)b\ = c 



(4) 



the adjoint operator effects the reversed process. The 
corresponding diagonal elements are build up by re- 
spectively . Now the new Hamiltonian can be written 
as: 



N-l 



gn{a) = a [n. 
Vn{a) = a 



A B 



(5) 



with the Hamiltonian H from Eq. 
rate for the CR process. 



), and a labels the 
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III. CALCULATIONS AND RESULTS 



In order to investigate the dependence of the viscosity 
on the CR rate a we performed DMRG calculations for 
the rates 1/iV^ and The choice of these rates 

is based on two estimates. The first is very simple: A con- 
straint is released when an end repton retracts back into 
the tube. The probability that the constraint is exerted 
by an end repton is 2/iV, provided that the probability 
for the presence of a segment is equally distributed. The 
hopping of a single repton is a process of rate one and so 
we find that the rate for a CR is of the order For 
the second estimate we assume that a CR is caused by a 
tube renewal. Later in this paper we will show the cal- 
culation of this estimate within the Rouse model which 
leads to the result that the rate a is proportional to the 
inverse relaxation time and scales therefore with \/N^. 
The viscosity which was calculated with these rates for 
the mechanism of AB exchange is plotted in Fig. ^. It 
can bee seen that the viscosity decreases with increasing 
CR rate which is reasonable because the relaxation of the 
tube is accelerated by this process. For the investigation 
of the scaling it is more useful to plot the local slope or 
the effective exponent 



ZN 



ln(Af-h 1) -ln(iV- 1)' 



(6) 



against 1/\/N as introduced in Ref. ^. The jdioice of the 
abscissa is motivated by the formula of DoiO which pre- 
dicts a correction to the TV ■^-scaling in the order of 1/ \/N 
by taking tube length fluctuations into account. It should 
be mentioned that in this paper N labels the number of 
segments which is proportional to the length of the poly- 
mer while in Ref. ^ it labels the number of reptons which 
is the number of segments plus one. Surprisingly the lat- 
ter interpretation shows a better agreement with the Doi 
formula while the former shows a better agreement with 
the experiments since the range where zn ~ 3.3 ± 0.1 is 
much broader. 

The influence on the effective exponent zn is plotted in 
Fig. Ij. The influence is non-monotonic in a: While 
the rate 1/A'^ causes an obvious shift down in compari- 
son with the data without CR, the data of rate 1/A^^ is 
located above the curve without CR. The data of rate 
1 /N^ is also located above the curve without CR but at 
a smaller distance. 

We remark that from a theoretical point of view also 
the rate a = 1 is interesting. If A- and i?-particles ex- 
change with the same rate as particles with vacancies, the 
particle can diffuse freely, as they feel no longer restric- 
tions. This means that the reptation model transforms 
into the Rouse model as the tube can change its form 
freely. This is why rate a = 1 represents a possibility 
to verify the model: The relaxation time should scale in 
the limit A^ ^ oo as A^^. But the DMRG calculation 
of the effective exponent with this rate yields a diverging 
curve. On the other hand, if the second excitation is con- 
sidered, one can see that it scales as A~^. This suggests 




FIG. 4: The viscosity [r] oc E]^^) for several CR rates of AB 
permutation. The solid line shows the result for the standard 
RD model without CR. 



without CR 
' Rate - 
' Rate - l/N" 
' Rate ~ l/N 




FIG. 5: The effective exponent zn for several CR rates of AB 
permutation. The solid line shows the result for the standard 
RD model without CR. 



that the first excitation is caused by an other relaxation 
mechanism and the second excitation yields the desired 
time, which scales with the expected exponent whereby 
the transition to Rouse dynamics is verified. 
The bounds 1/A > a > are rather weak and a self- 

consistent approach is needed for determining the appro- 
priate scaling of the rate a. As a CR event is caused 
by the tube renewal of a polymer, this process is not in- 
dependent of the relaxation time. The relaxation time 
itself is affected by the CR rate. So the rate is physical, 
if the relaxation time, which is calculated using this rate, 
yields the same CR rate. 

To calculate the self consistent rate asc we proceed as 
follows: The dependence of the CR rate on the relaxation 
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t(a) 



without C.R. 
S.C. 

DMRG with a,^ 



FIG. 6: Relaxation time versus CR rate a; asc is given by 
the intersection point of the curves rjv and r (here exemplified 
for = 30). 



FIG. 7: The new DMRG calculation of viscosity - a good 
coincidence with the data of the intersection points can be 
seen. 



time air) is estimated by an analytical calculation using 
Rouse dynamics. The DMRG calculation with the rates 
mentioned above yields directly the dependence of the 
relaxation time on the CR rate rjv(a) for the respective 
chain lengths. The inverse function of q;(t), T{a), is com- 
pared with TN^a): The intersection of the curves yields 
the self consistent rate asc in first approximation. 
First we calculate the dependence of the CR rate on the 
relaxation time. The polymer whose relaxation time we 
want to calculate finally is hindered in his free move- 
ment by an other polymer. The position which exerts 
the constraint may have the distance s from one end of 
the constraining polymer. So, the constraint is released 
when the polymer moves either the distance s into the 
one direction or the distance N — s into the other direc- 
tion, where N labels the length of the polymer. In Ref. 
^ an expression is presented which provides the distribu- 
tion i'{s,t) of the probability that the tube segment s 
was at time t not yet reached by the ends: 

p.odd 

with the relaxation time r. The calculation of this for- 
mula bases on the assumption that the polymer diffuses 
in between the tube due to Rouse dynamics. The comple- 
mentary probability distribution (/)(s,t), which indicates 
the probability that at time t the segment s is already 
reached by one of the ends is then 

4>{s,t)^\-ij{s,t). (8) 

The derivate with respect to time 




is the first passage time density, i.e. the probability per 
time that the segment s is reached by one end just at 
time t. Now we can specify the mean first passage time 

POO 

fi{s)=. / dttf{s,t) (10) 
Jo 

which indicates how long it takes in average until segment 
s is reached by an end. Finally we average over all s 
provided that each segment builds up an entanglement 
with the same probability 

1 /-^ 

= ^ 7o rfsAi(s)- (11) 

This is the desired time: In average the release of a con- 
straint will take the time the CR process takes place 
with rate a = ■ The calculation of the integrals is 
elementary mathematics and we only show the result: 

^ = 0.822r a = 1.22t"^ (12) 

This calculation is based on a continuous description 
of the polymer. However, we are considering a lattice 
model of reptation so that boundary effects might 
not be taken into account correctly by the continuum 
expression (^l]). Therefore we performed the above 
calculation as well for the discrete case. As the result 
for the discrete case converges quickly to the continuous 
case but consists of only numerical evaluable series, we 
restricted ourselves to the continuous case here. There is 
an error of less than 2 per cent for as little as 10 reptons. 

Fig. ^ shows the curve T{a) — 1.22a^^ and the re- 
laxation times for the CR rates a — and 
exemplified for N = 30, the data point for a = 1/N 
is out of the plotted range. In order to investigate the 
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FIG. 8: The effective exponent zn of the viscosity with seff- 
consistent CR. 



dependence of the relaxation time on the CR rate we 
interpolated the data points linearly by fitting them 
with splines. The data points are is well approximated 
by this interpolation in the region where the intersection 
point is expected. In this way we determined the self 
consistent relaxation times and rates for the lengths 
TV = 8 ... 50. As the lengths = 4 and TV = 6 do 
not lead to reasonable intersection points we did not 
take them into account. To verify the self consistence a 



new DMRG run was done with a fit of the determined 
rates; as can be seen in Fig. Q the resulting relaxation 
times coincide well with the ones determined by the self 
consistence condition, so that the linear interpolation 
represents a sufficiently good approximation. 
In Fig. H the effective exponent is shown for the self con- 
sistent relaxation times and for the DMRG calculation. 
The exponent is shifted to higher values in the region of 
small N ^ but a broadening of the region where ~ 3.3 
could not be observed. 



The Rubinstein model is not microscopic and hence 
there is some freedom in implementing rtiicroscopic pro- 
cess such as CR. A different mechanismll3 leads to the 
creation and annihilation of particles in the bulk. We 
performed similar analysis for this mechanism with qual- 
itatively similar results.Ej 

We conclude that Rouse-based calculations do not lead 
to a proper description of crossover behavior in the vis- 
cosity of polymer melts when Rouse theory is used as 
a self consistency input in the mesoscopic and generally 
quite successful Rubinstein model for reptation. We can- 
not rule out that a fully self-consistent implementation 
of CR (without Rouse assumption) leads to a crossover 
regime closer to empirical evidence, but the broad range 
of CR rates studied here suggests that bulk shape fluc- 
tuations of the tube which result from constraint release 
do not significantly broaden the crossover range of the 
viscosity. 
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